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SUMMARY 

The  maximum  signal -to-noise  ratio  array  processor  is  formulated 
in  beam  space.  Expressions  for  the  optimum  narrowband  weight 
vectors  and  array  gain  are  derived.  Some  general  properties 
of  the  beamspace  formulation  are  derived  and  conditions  for  the 
equivalence  of  the  array  and  beam  space  formulations  are  proved. 

Examples  using  both  simulated  and  sonar  data  are  given  to  compare 
the  beam  and  array  space  formulations. 
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1.  INTRODUCTION 

The  performance  of  an  array  bcamformer  may  often  be  improved  by  multiplying  the 
receiver  outputs  by  a  set  of  weighting  (or  shading)  coefficients  prior  to  addition. 
In  particular  Edelblute,  Fisk  and  Kinnison (ref . 1)  have  derived  a  set  of  weights 
which  for  a  known  noise  field,  maximise  the  output  signal -to-noise  ratio  of  a 
beam  steered  in  a  given  direction.  These  weights  are  determined  by  the  crosspower 
spectral  matrix  of  the  receiver  noise  outputs  and  hence  reflect  the  angular 
distribution  of  noise  sources. 

In  practice  the  receiver  outputs  may  not  be  accessible  and  so  an  estimate  of  the 
receiver  noise  CTOsspower  spectral  matrix  cannot  be  directly  obtained.  One 
important  practical  case  is  where  only  beamformed  outputs  (ie  conventional  beams) 
are  available.  A  set  of  weighting  coefficients  which  maximise  the  output  signal - 
to-noise  ratio  are  derived  in  this  paper  from  conventional  beam  rather  than  receiver 
outputs.  Expressions  for  the  array  gains  are  also  derived.  The  expressions 
derived  are  a  generalisation  of  the  equations  obtained  by  Vural(ref.2)  using  a  point 
constrained  minimum  power  criterion. 

In  general  the  array  space  and  beam  space  formulations  are  not  equivalent.  However 
as  the  number  of  independent  beams  used  in  the  beam  space  formulation  increases  the 
beam  space  gains  approach  the  array  space  gains.  In  Section  3  the  equivalence  of 
the  two  approaches  is  proved  under  some  general  conditions. 

Some  general  properties  of  the  beam  space  formulation  are  derived  in  Section  4. 

In  Section  5  some  theoretical  and  simulated  examples  are  presented  to  illustrate 
the  relationship  of  the  conventional ,  the  optimum  array  space  and  the  optimum  beam 
space  processors.  The  simulation  examples  illustrate  some  practical  limitations 
on  the  gains  that  can  be  achieved  using  the  beam  space  methods.  Finally  in 
Section  5  some  examples  of  the  application  of  the  beam  space  optimisation  to  sonar 
data  are  given. 

This  work  is  part  of  a  continuing  R$D  programme  in  signal  processing  for  underwater 
detection  and  has  been  carried  out  under  task  DST  79/069. 


2.  DERIVATION  OF  OPTIMUM  BEAM  WEIGHTS 

2.1  Notation 

Let  y  denote  the  vector  of  M  conventional  beam  outputs  at  a  frequency  f,  and 
let  y  be  a  weighted  scalar  sum  of  these  beam  outputs,  ie 

H 

y  =  w  y 

where  w  is  the  (in  general  complex)  vector  of  weights. t  If 


y  =  s  +  n 

where  and  (i=l,2,...M)  are  the  signal  and  noise  components  in  the  i-th 

beam,  then  the  crosspower  spectral  matrices  of  the  signal  and  noise  beam 
components  are  defined  by 

(0) .  .  =  <  s . s .  > 

ij  i  J 


t  The  superscript  denotes  the  Hermitian  transpose  of  a  matrix  or  a  vector. 
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and 


<  n.n.  > 
i  3 


respectively,  where  <  >  denotes  ensemble  averaging  and  i,j  =  1,2,...,M  where 
M  is  the  number  of  (preformed)  beams  used. 

If  the  signal  is  a  narrowband  plane  wave  (of  unity  amplitude)  then 


Q 


=  ss 


il 


where 


S  =  VH  V(0), 
(V)kj  -  e-27lifrkj 


and 


vk(») 


e-2.ifrk(0) 


where  k=l,2,...,K  and  K  is  the  number  of  receivers. 


The  time  delay  t 


kj 


is  the  time  delay  relative  to  some  arbitrary  reference 


point,  at  the  k-th  receiver  of  a  plane  wave  from  th  j-th  direction  (corres¬ 
ponding  to  one  of  the  chosen  steer  directions  for  the  M  preformed  beams).  The 
delays  t^(0)  are  the  time  delays  corresponding  to  the  steering  direction,  0} 
of  the  optimally  weighted  beam. 


The  signal -to-noise  ratio  in  the  derived  beam,  y,  is  defined  as 


SNR 


H 

w  Qw 
IL 

w  Rw 


Array  gain  is  defined  by 


(SNR)  beam 

(SNR)  omnidirectional  receivers. 


2.2  Derivation 


It  can  readily  be  shown  that  the  SNR  in  the  derived  beam,  y, 
choosing  w  such  that 


IL 

w  Qw 
IL 

w  Rw 


Rw 


is  maximised  by 


Qw 


to 
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Equation  (1)  is  satisfied  if 


XfT  1  s 


where  X  is  any  constant.  imposing  the  constraint  of  a  fixed  (ie  unity)  response 
in  the  steering  direction  defined  by  v ^  (0 )  implies  that 


w  s  =  1  . 

o 


This  constraint  is  imposed  to  ensure  that  the  signal  power  in  the  optimum  beam 
increases  as  the  signal  power  increases.  Consequently  equation  (2)  becomes 


5%  1  s 


IT  1  VHv 

v1  W  1  VHv  .  (3) 

As  discussed  in  Appendix  A  the  case  where  R  is  singular  and  V  is  of  full  rank 
can  also  be  treated  and  an  appropriate  generalisation  of  equation  (3)  is  to 
replace  R" 1  by  R  ,  the  Moore  Penrose  pseudoinverse (ref . 3)  which,  in  addition 
to  maximising  the  output  SNR  also  minimises  the  superdirectivity  of  the  weights. 

It  directly  follows  from  the  generalised  version  of  equation  (3)  that  the 
maximum  output  SNR  for  a  unit  amplitude  signal  is  given  by 

(SNR)  =  sHR+s  . 


For  these  optimally  weighted  beams  the  expression  for  array  gain  reduces  to 


H  + 

y  _  sRs 


where  Nq  is  the  omnidirectional  noise  power  at  any  receiver.  (The  assumption 
of  a  homogeneous  medium  and  identical  receivers,  ie  no  shading,  has  been  made). 

Furthermore  if  R,  the  beam  noise  crosspower  spectral  matrix  has  been  derived 
from  normalised  receiver  outputs  ie,  the  noise  output  of  the  j^*1  receiver, 

(x) 

n  ,  is  transformed 
J 


J  J 


'  )  J 
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then  the  optimum  array  gain  in  the  steer  direction  specified  by  s  is 


H  + 

s  R  s. 


3.  SOME  EQUIVALENCES  OF  THE  ARRAY  SPACE  AND  BEAM  SPACE  FORMULATIONS 

lor  narrowband  signals  the  conventional  beam  outputs  y^ ,  i=l,2,...M  can  be 
related  to  the  receiver  outputs  x^,  i=l,...,K  by  the  equation 


i  r*  * 

y  =  V  x 


(4) 


where  V  is  defined  in  the  previous  section.  The  crosspower  spectral  matrix  of 
the  M  beam  outputs,  R^,  is  defined  by 


(R  (f)n^y)*  (f)  > 

y  ij  i  j 


for  i , j =1 , 2 , . . . ,M  and  the  crosspower  spectral  matrix  of  the  K  receivers  is 
defined  by 


(R  ) .  .  =  <  nfx)  (f)n^ * (f)  > 

X  lj  1  J 

for  i , j=l ,2, . . . ,K  where  nfX^  and  are  the  i-th  normal isedt  receiver  noise 

outputs  and  j-th  beam  noise  outputs  respectively. 

Assuming  the  signal  and  noise  are  uncorrelated  it  follows  from  equation  (4)  that 


R  =  V  R  V. 
y  x 

From  reference  1  the  optimum  weight  vector,  wq,  using  receiver  outputs  is 


IT  *  v 
x 

vAr'v 

x 


(5) 


and  the  corresponding  array  gain,  g^,  is  given  by 


X  ri  r  , 

g  =  v  IT  v  . 

60  X 


(6) 


i 


j 


+  The  use  of  normalised  receiver  outputs  as  defined  in  the  previous  section  is 
simply  a  convenience  which  allows  simple  expressions  for  the  crosspower  spectral 
matrices  and  a  concise  formulation  of  array  gain  within  beam  space. 


5 
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y 

From  the  previous  section  the  optimum  weight  vector,  w' ,  using  beam  outputs  is 


R+s 

y 

n* 

s  R  s 

y 


(7) 


where  s  = 


V  v. 


Also  the  corresponding  array  gain,  g^,  is  given  by 


g 


y 

o 


R 

y 


(8) 


Two  conditions  under  which  equations  (5)  and  (7)  and  (6)  and  (8)  become  identical 
will  now  be  proved. 

Case  A:  K  =  M  and  V  nonsingular. 

Since  is  assumed  nonsingular  it  follows  that 


It  then  follows  that 


r+  =  ir 1  =  v  '  ir  1  v  1 H 
y  y  * 


sVs 

y 


v  v  \r  1  i r  1  v  '  v  v 

X 


-  vV  V 
X 


Also  the  optimally  weighted  beam  output  is  identical  to  the  optimally  weighted 
array  output,  ie 


yH 
w'  y 
o  7 


H  ..  +H 
v  VR  y 

y 

v^V  R+  VHv 

y 


V«v  V  1  IT  1  V" 1  fl  VH  x 
x 

vHy  |  I  y"  1  H  yFl 
X 


WSUI.-01 7i)  -TR 


0 


H  _  , 

v  ir  x 


v  ir  1  v 

X 


xH 

=  w  x . 
o 


Ihus  for  K  «  M  and  V  nonsingular  the  two  techniques  are  equivalent, 
base  B:  M  >  K  and  V  is  of  rank  K 

In  this  case  only  k  of  the  M  beams  formed  arc  independent.  For  most  arrays  this 
can  be  satisfied  by  choosing  the  M  steering  directions  to  be  distinct (ref .4) . 
Since  V  has  K  independent  rows  it  fol lows (ref . 3)  that 


v+  =  vH(wHr 1 


and  consequently 


VV 


=  v+HvH  =  I 


\s  a  result  it  can  readily  be  proved  that 


R+  =  V+  If  1  V 


+H 


it  then  follows  that 


11  D  + 

s  R  s 

y 


=  V 


’V  V+  (T  1  V+H  VH 


11  1 
v  IT  v 

X 


In  a  similar  manner  to  the  previous  section  the  optimally  weighted  receiver  and 
beam  outputs  can  also  be  shown  to  be  identical,  ic 

y\\  1  H  +H 

w/  y  =  - v  V  R  y 

o  7  y  y 

ir. 


=  1  H  D-  ,  .  +H  H 

-  v  V  V  R  V  V  x 

X 

*0 
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g 


o 


x 


These  examples  imply  that  if  the  number  of  independent  beams  is  greater  than  or 
equal  to  the  number  of  receivers  then  in  principle  it  is  possible  to  achieve  the 
full  array  space  gains  using  the  beam  space  approach.  Furthermore  as  discussed 
in  reference  4  choosing  the  steering  directions  to  be  distinct  implies  linear 
independence  of  the  steering  vectors.  This  implies  that  at  any  frequency  the 
full  array  gains  can  be  achieved  by  increasing  the  number  of  beams  used.  In 
practice  the  finite  precision  used  in  the  calculations  prevents  this  and  in 
Section  S  some  examples  are  given  to  illustrate  this  difficulty. 


4.  PROPERTIES 

In  this  section  some  general  properties  of  the  beam  space  formulation  will  be 
detailed . 

4.1  Gains 

y 

The  array  gains  using  the  beam  space  formulation  ie  the  g'  are  always  less  than 
or  equal  to  those  derived  from  the  array  space  formulation  ie  gQ . 

The  proof  follows  from  Rao(p.48  ref. 5): 


x  H  -  !  H  H 

p  s  v  IT  v  =  max  x  vv  x 

°  X  xe  F.  H 

K  x  R  x 
x 


where  is  the  vector  space  of  dimension  K.  Thus  for  any  vector  ye  it 
follows  that 


H  H 

y  vv  y  <  0X 

6o 

yHR  y 

J  X  / 


In  particular  take 


y  =  VR+  v"v 

y 


it  then  follows  that 
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From  the  definition  of  R  and  the  Moore-Penrose  pseudoinverse  it  follows  that 

y  1 1  I  {  ^ 

the  I. US  of  the  inequality  reduces  to  v  VR  V  v  and  hence  v?  ^  g  X.  The 

y 

conditions  for  the  equality  to  hold  have  been  investigated  in  the  previous 
sect  ion . 


•1.2  iinin  as  a  function  of  number  of  input  beams 

In  a  previous  classified  report  the  author  has  shown  that  in  array  space  the 
optimum  array  gain  estimates  always  increase  with  the  addition  of  extra 
receivers.  This  is  not  the  case  in  conventional  beamforming  where  at  low 
frequencies  the  addition  of  extra  receivers  may,  provided  the  array  aperture 
is  not  increased,  actually  result  in  a  poorer  array  gain. 

\  similar  result  proved  in  Appendix  II  and  generalising  the  previous  derivation, 
holds  for  the  beam  space  formulation  ie,  optimum  beam  space  gains  in  general 
increase  when  additional  input  beams  are  used.  However  the  array  gains  attain 

a  r.uxinum  when  the  number  of  beams  used  equals  the  number  of  receivers  and  all 
these  beams  are  independent.  As  shown  in  Appendix  II  when  R  formed  from  using 

iv  1  beams  becomes  singular  (or  when  the  rank  of  R  is  not  increased  by  the 

additional  beam)  then  the  optimum  beam  space  gain  is  no  longer  increased  by  the 
addition  of  an  extra  beam. 


5.  EXAMPLES 


In  thv  section,  array  gain  estimates  for  a  conventional  beamformer  will  be  denoted 
is  g  and  can  readily  be  proved  to  be  given*  by 


8C 


K 


2 


vHR  v 

X 


5.1  Uncorrcl at ed  receiver  noise 

In  this  case  the  receiver  noise  covariance  matrix  is  defined  by 


R  »  I . 
x 


It  trivially  follows  that 


g 


c 


K 


for  all  directions. 

Furthermore  e?  reduces  to 
o 


H  ♦  H 
Y(V  V)  V  v 


*  I  he  assumption  that  the  signal  lies  in  the  look  direction  is  made  henceforth. 


I 


-  9  - 
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and  consequently  the  optimum  beam  space  gains  will,  as  would  he  expected, 
depend  on  the  number  of  beams  used  and  the  steer  direction. 

(a )  Single  beam 

Using  just  a  single  beam  steered  in  a  direction  corresponding  to  the 
i-th  column  of  V  then 

R  -  vV  v.  =  K 
y  i  l 

where  v.  is  the  i-th  column  of  V. 

i 

Thus 


T 

..  j 


gy0  =  v(0)  v(0) 

which  is  just  the  polar  diagram  of  the  array  when  steered  in  the  i-th 
direction.  Since  only  one  beam  has  been  used  to  effect  the  optimisa¬ 
tion  this  result  is  hardly  unexpected. 

(b)  M  orthogonal  beams  (M  <  K) 

In  this  case 


v'V  .  K  IH  .  >  (v"vr '  ■  i  IM 

and  hence 

8o  =  £  vH(0)  wV*) 

M 

1  H...  H  ... 

=  F  )  v  (0)  v.v  v (0) 

K  L,  J  J 


ie  the  sum  of  the  M  polar  diagrams. 
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Thus 


10 


H 


y 

o 


v'V 


1 

0 


0 


=  K. 

Thus  the  optimum  beam  space  gains  in  the  M  directions  determined  by  the 
v^fs,  i  =  1,  2,  ...  M  will  be  equal  to  the  optimum  array  space  and 

conventional  array  gains. 

5.2  Simulation  results 

To  investigate  the  dependancc  of  the  beam  space  gains  on  the  number  of  input 
beams  used  the  crosspower  spectral  matrix  was  simulated  for  a  number  of  noise 
fields  and  the  gains  were  evaluated.  The  conventional  array  gain  and  the 
optimum  array  space  gain  were  also  evaluated.  An  array  of  two  concentric 
rings  each  of  5  receivers  (.sec  figure  1)  was  used  for  the  simulation. 

(a)  The  array  gains  for  a  unit  receiver  noise  covariance  matrix  (corres¬ 
ponding  to  uncorrelated  receiver  noise)  are  plotted  in  figure  2  where 
r/X  =  1/8.  The  number  of  beams  used  in  deriving  the  beamspace  gains 

was  varied  from  1  to  5;  in  each  case  the  beams  were  equi spaced  from 
0  to  2  7 r.  For  M=1  the  polar  diagram  of  the  array  is  apparent.  For 

M  >  1  a  gain  of  10  log  10  (sec  previous  section)  is  obtained  whenever 
the  steering  direction  coincides  with  one  of  the  preformed  beam 
directions.  Furthermore  at  this  frequency  only  5  beams  are  needed 
to  achieve  the  full  array  gain  of  10  log  10. 

The  limiting  number  of  beams  needed  is  a  function  of  frequency  and  are 
plotted  in  figure  3.  This  limiting  number  was  calculated  as  the 
minimum  number  of  beams  used  such  that  the  beamspace  gains  lay  within 
an  average*  of  h  dR  of  the  limiting  array  space  gains.  For  such  a 
noise  field  this  gives  the  maximum  number  of  beams  that  need  be  used 
at  any  given  frequency. 

I b)  Uncorrelated  receiver  noise  and  a  single  interference 

The  array  gains  for  the  case  where  the  above  noise  field  has  also  a 
10  dB  interference  at  180  are  plotted  in  figure  4  where  the  analysis 
parameters  are  the  same  as  figure  2.  The  number  of  beams  necessary 
to  achieve  the  limiting  array  gains  are  plotted  as  a  function  of 
frequency  in  figure  3.  Comparing  the  two  plots  of  figure  3  it  can  be 
seen  that  the  interference  requires  roughly  an  additional  beam  (particu¬ 
larly  at  very  low  values  of  r/X)  to  be  used  in  order  to  achieve  the  full 
array  gains. 

(c)  Isotropic  noise 

for  an  isotropic  noise  field  the  cross -correlation  between  any  two 
receivers  separated  by  a  distance  of  d  is  given  by 

sin  2*d/X 
2*d/X 


averaged  over  the  azimuth  angles  0-2* 
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The  array  gains  for  the  array  of  figure  1  in  such  a  noise  field  are 
plotted  in  figure  5  where  the  number  of  beams  used  was  varied  from 
1  to  5  and  r/X  =  1/8.  The  number  of  beams  necessary  to  achieve 

the  limiting  array  gains  are  plotted  as  a  function  of  r/X  in  figure  3. 
Apart  from  very  low  values  of  r/X  the  values  are  the  same  as  for  the 
uncorrelated  receiver  noise  case. 

(d)  Isotropic  noise  and  a  single  interference 

As  discussed  earlier  it  is  in  principle  always  possible  to  obtain  the 
full  array  gains  by  increasing  the  number  of  beams  used  in  the 
bcamspace  approach.  However  in  this  ease  it  was  found  that  as  the1 
number  of  beams  increased  the  factor i sat iont  of  the  beam  crosspower 
spectral  matrix  failed. 

A  method  that  was  partially  successful  in  overcoming  this  difficulty 
was  to  perturb  the  diagonal  elements  of  Ry  by  a  small  positive  number 
(1(T5).  This  permitted  factorisation  when  r/X  >  1/6  and  in  these 

cases  the  full  array  gains  were  achieved  using  either  9  or  10  beams 
depending  on  frequency.  Increasing  the  perturbation  ( 1  (T  4 )  permitted 
factorisation  at  most  lower  values  of  r/X.  However  in  general,  the 
limiting  beamspace  gains  (at  r/X  =  1/8)  were  significantly  lower  than 

those  achieved  in  the  array  space  and  often  as  the  example  of  figure  6 
shows  the  limits  could  be  achieved  using  less  than  10  beams  (eg  5  in 
figure  6) . 

The  inference  from  these  results  is  that  all  10  beams  are  needed  at 
all  frequencies  to  achieve  the  maximum  gains.  At  lower  frequencies 
the  ill-conditioning  of  the  matrices  limits  the  gains  that  can  be 
achieved  in  practice.  However  if  uncorrelated  receiver  noise  is  also 
present  (and  in  practice  this  certainly  will  always  he  the  case)  then 
the  limiting  gains  may  be  achievable.  If  -10  dB  (re  the  isotropic 
noise  power)  of  uncorrelated  receiver  noise  is  added  in  the  above 
example  then  figure  7  is  obtained  and  it  can  be  seen  that  using  10 
beams  allows  the  full  array  gains  to  be  achieved. 

5.3  Application  to  real  data 

The  methods  described  have  been  used  to  estimate  the  optimum  gains  from  the 
array  described  in  the  previous  examples  when  processing  sonar  data.  A 
block  diagram  of  the  processing  is  given  in  figure  8.  After  narrowband 
filtering  the  data  (by  means  of  the  FFT)  to  a  5  Hz  binwidth,  the  receiver 
crosspower  spectral  matrix,  R^ ,  was  estimated  and  normalised  at  each  frequency 

of  interest.  This  was  factored  and  the  factored  matrix  was  used  to  estimate 
the  conventional  and  optimum  array  space  gains.  The  factored  matrix  was  also 
premultiplied  by  the  matrix  (V)  of  phase  delays  whose  columns  corresponded  to 
the  steering  directions.  From  this  R  was  derived  (see  figure  8),  factored 

y  y 

and  inverted  and  hence  the  array  gain  estimates  g^  were  evaluated. 

The  array  gain  estimates  as  a  function  of  azimuthal  angle  for  the  array  of 
figure  1  are  shown  in  figures  9  and  10  for  r/X  equal  0.125  and  0.42  respectively. 
Since  one  of  the  outer  ring  hydrophones  was  faulty  only  9  hydrophones  were  used 
in  the  analysis.  The  data  also  contained  random  phase  errors  and  this  tended 
to  bias  the  optimum  array  gain  estimates  positively.  In  both  these  examples 
no  problems  were  encountered  in  the  beam  crosspower  matrix  factorisation 
although  at  other  lower  values  of  r/X  and  also  when  using  arrays  with  a  larger 
number  of  hydrophones  on  the  same  data  the  matrix  factorisation  failed  and  the 
perturbation  technique  of  the  previous  section  had  to  be  used. 


t  The  matrix  factorisation  arithmetic  was  effected  in  double  precision  on 
single  precision  inputs  on  a  4  byte  word  machine. 
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The  two  examples  were  chosen  to  illustrate  two  different  environments;  in 
the  first  example  (.figure  9)  the  noise  field  exhibited  only  a  weak  azimuthal 
dependance  whereas  in  figure  10  the  presence  of  a  strong  interference  at 
264  implied  a  strong  azimuthal  dependance  of  the  noise  field.  Both  results 
showed  an  improvement  in  the  beam space  array  gain  estimates  as  the  number  of 
beams  used  was  increased.  When  the  number  of  beams  equalled  the  number  of 
receivers  (figure  9)  then  the  full  optimum  array  space  gains  were  achieved  in 
both  cases.  Both  examples  showed  distinct  dips  in  the  gain  estimates  when 
only  a  small  number  of  beams  were  used  (eg  3  in  figure  9  and  4  in  figure  10). 

Care  should  be  taken  to  avoid  associating  these  dips  with  an  azimuthal 
variation  of  the  noise  field;  they  are  simply  a  result  of  using  an  insufficient 
number  of  beams.  Figure  10  differs  from  figure  9  in  that  significant  improve¬ 
ment  over  the  conventional  gains  can  be  achieved  by  using  only  a  few  beams, 
eg  4,  whereas  in  figure  9,  the  gains  resulting  from  using  just  3  beams  arc 
often  significantly  lower  than  the  conventional  array  gains.  This  is  because 
only  a  single  beam  pointed  ac  or  near  the  interfering  source  is  sufficient 
to  partially  null  it  and  hence  produce  the  significant  array  gains  of  figure  10. 

Two  further  features  are  of  particular  interest  in  figure  10.  The  first  is 
the  depth  of  the  null  as  a  function  of  the  number  of  beams  used;  this  depth 
is  a  good  indication  of  the  sensitivity  of  the  technique  to  receiver  errors. 

As  the  number  of  beams  increases  the  depth  of  therull  decreases  indicating 
increased  sensitivity  as  the  number  of  beams  used  increases.  The  second 
feature  also  of  interest  is  the  distinct  dip  at  68  in  the  array  gain  estimate 
for  the  6  beam  case.  Tt  is  common  to  associate  such  dips  in  array  gain 
estimates  with  interfering  targets,  however  as  indicated  by  both  the  conventional 
and  optimum  array  space  gains  there  is  no  interference  at  68  .  Thus  as 
discussed  above  for  the  three  and  four  beam  case  this  dip  must  be  a  result  of 
using  an  insufficient  number  of  beams.  This  again  illustrates  the  importance 
of  choosing  the  number  of  beams  correctly  if  the  array  gain  estimates  are  to 
be  used  for  the  detection  of  interferences. 


<>.  SUMMARY 

I  he  optimum  signal -to -noise  ratio  processor  has  been  formulated  using  only  narrowband 
preformed  beam  outputs.  The  relationship  of  this  to  the  array  space  formulation 
has  been  investigated,  and  in  priciple  the  full  array  space  gains  can  be  recovered 
from  the  beamspacc  formulation  provided  the  number  of  beams  is  chosen  appropriately. 

A  number  of  theoretical  properties  together  with  examples  illustrating  these 
properties  have  been  given.  Some  practical  limitations  of  the  beamspace  methods 
are  also  illustrated.  The  method  has  been  applied  to  sonar  data  from  an  array  of 
hydrophones  to  verify  its  viability  and  demonstrate  some  of  the  theoretical 
predictions.  The  results  indicate  that  a  trade-off  of  increased  sensitivity 
versus  more  uniform  (in  an  angular  sense)  array  gain  estimates  must  be  made.  This 
problem  of  choosing  the  correct  number  of  beams  is  receiving  further  attention. 
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APPENDIX  I 

DERIVATION  OF  MAXIMUM  SIGNAL-TO-NOISE  RATIO 
PROCESSOR  FOR  A  SINGULAR  BEAM  COVARIANCE  MATRIX 

Let  he  a  positive  definite  matrix  and  define  R^  by 

Ry  =  V1'  KxV  (1.1) 

where  V  is  an  K  x  M  matrix  and  M  >  K. 

Since  R^  is  a  nonsingular  it  follows  that 


H 

u 


R  u 
x 


0  =  >  u 


0  . 


(1.2) 


However  from  equation 
u'  such  that 


(1 . 1)  the  rank  of  R 


y 


is  K  and  hence  there  exist 


(M-K)  vectors 


•H  ? 

u  R  u  =  0  . 

y 

Substituting  equation  (1.1)  in  the  above  expression  it  follows  from  equation  (1.2) 
that 


Vu f  =  0  . 

The  general  solution  to  the  above  equation  is(rcf.3) 

u'  =  (1  -  V  +  V)  7. 


where  z  is  arbitrary  and  V*  is  the  Moore-Penroso  pseudoinverse  of  V. 

Let  s  (0)  =  V%(0)  be  the  steering  vector,  it  thus  follows  that 

SH  u'  =  vH(0)  V(I-V+V)z 

=  0 

since  VV+V  =  I  . 

Thus  it  follows  that  any  vector  which  completely  nulls  out  the  noise  power  also 
nulls  out  the  signal  power. 

This  important  conclusion  now  allows  a  derivation  of  the  maximum  SNR  processor  for 
Ry  singular  to  be  effected.  That  is  the  expression  for  the  SNR  in  a  given  output 

beam  ie 
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II  I  i 

w  ss  w 

(SNK).  -  — - 

'beam  I L 

w  U  w 

y 

can  be  maximised  subject  to  the  constraint  of  a  finite  response  (say  unity)  in 
the  look  direction. 

This  is  a  formulated  equivalently  as  minimising  R^w  subject  to  the  constraint 

sHw  =1  (1.3) 

Solving  in  the  usual  way  by  introducing  a  Lagrangian  multiplier  and  differentiating 


with  respect  to  w  implies  that 


R  w  =  Xs  . 

y 


Any  solution  of  this  equation  can  be  expressed  in  the  form 

w  =  XR+s  +  (I  -  R+  R  )z 

y  y  y 

where  z  is  arbitrary  and  R^  is  the  pseudoinverse  of  R^. 

lor  the  case  of  M  >  K  but  V  of  full  rank  it  can  be  shown  that  (see  Section  3) 


R  =  V  fT  1 V 

y  x 


and  consequently  that 


Thus  w  =  XR+s+(l-VV)z  and  so 

y 


AsHR+s  +  vflV(I  -  V+V) 

y 


=  AsVs  . 


Choosing  X  such  that  equation  (1.3)  holds  then  implies  that  a  possible  solution  for 


w  is 
o 


and  consequently 


s"  R+s 


(SNR)  =  s“  R+s  . 

v  J  max 


WSkl.  0197-1 K 


16 


APPENDIX  II 

INCREASE  IN  OPTIMUM  ARRAY  CAINS  BY  Till!  AUDITION  OF  AN  EXTRA  INPUT 


Let  R^  denote  the  crosspower  spectral  matrix  formed  using  n  'normalised'  input 

channels.  These  could  be  either  array  elements  or  conventional  beams.  Fet 
R  j  be  the  crosspower  spectral  matrix  formed  using  the  same,  n,  input  channels 

plus  an  additional  channel.  It  follows  that  R  ,  can  be  written  in  the  form: 

1  n+1 


"„.i  ■  ( \b\ 

V  l>H  d/ 


where 


<  X  ,  X.  >  , 

n+1  i 


d 


n  +  1 


> 


and  the  X.'s  are  the  Fourier  transforms  of  the  'normalised'  i  channel  outputs. 
i 

If  u  is  any  'steering  vector'  for  n  input  channels  then  the  optimum  array  gain  is 
given  by 


where  R*  is  the  Moore-Penrose  pseudoinverse  of  R  .  Similarly  for  the  n+1  input 
n  n 

channels  the  array  gain  is 


g 


n+1 

o 


<“H  «>  c, 


(II. 1) 


where  a  is  the  complex  weight  required  by  the  additional  channel. 

Now  since  R  ,  is  positive  semidefinite  then 
n  +  1  1 


n+1 


=  A  A 


and  A  can  be  partitioned  A  =  (A  :  a)  where  A  is  n  x  n-1  and  a  is  n  x  1.  Thus 


n+1 


/  aHa  aH 

-  /  A  A  A  a 

Ha  II 
a  A  a  a 


I 


"1 
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It  can  easily  be  verified  that  theoremt  3.0.3  (p.00  of  Kao  and  Mitra)  a  1  so  holds 
for  the  pseudoinverse  and  hence 


where 

h  =  (A  A)  A  a 
=  A+a, 

and 

6  =  (a^a  -  AA*af  1  if  a  ^  M(A) 

=  0  if  af  M (A) 


where  M[A)  denotes  the  column  space  spanned  by  the  columns  of  A. 
Thus  it  follows  that  the  pseudoinverse  of  ^n+1  Is  given  by 


R  +  ,  =/  I k"  +  5R+  bbM  R+  '  6r+  b 

n+1  /  n  n  n  n 


-  6bH  R+ 
n 


Then  substituting  in  equation  (II.l)  it  can  be  shown  that 


gn+1  =  uH  R+u  +  6 (w-a*)  (w*-a) 


(II. 2) 


where  w  =  -u^  R*b .  Thus  provided  6  >  0  then  +  *  ^  u^R*u  anc*  s0  t^ie  arraX 
gain  is  never  decreased  by  the  addition  of  an  extra  input  channel. 

Now  since  K  ,  is  non-negative  definite  it  follows  that 

R  =  UH  a  ,  U 
n+l  iV  n+1 

where  IJ  is  a  unitary  matrix  and  A  ^  is  a  diagonal  matrix  such  that 


x.  =  (A  .) .  .  > 

i  v  n  +  l'n 


0. 


t  Two  typographical  errors  occur  on  p.66  of  reference  6,  a  is  an  m  vector  and 
a  e  M(A)  should  read  a  *  M(A). 


? 
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It  holds  that 


UH  A 


n+1 


U 


where  A  +  , 
n+1 


is  diagonal  and 


C  A 


n+1 ^  i  i 


>>  " 

1 

0  i  f  x  .  =  0 . 

i 


Thus  R  +  ,  is  non-negative  definite  and  so 
n+1 

xH  R+  t  x  >  0  V  x. 
n+1 

Taking  =  (0,  0,  ...»  1)  it  follows  that 

^n+Pn+1,  n+1  >  0, 

ie  5  >  0, 

The  ease  of  5  =  0  corresponds  to  the  case  of  where  the  rank  of  equals  that  of 

R  .  In  this  case  the  array  gain  is  not  increased  because  no  additional  informa¬ 
tion  is  contained  in  R  ,  that  was  contained  in  R  . 

n  +  1  n 

If  the  rank  of  R  ,  is  greater  than  that  of  R  then  5  >  0  and  hence,  in  general, 
the  array  gain  will  increase  since  w  V  a*  for  all  steering  vectors  u. 


max 


T 

t 

I 
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Figure  3  ; 


Figure  3. 


Limiting  number  of  beams  as  a  function  of  r/X 
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—  Conventional 

□  Optimum  array  space 
O  Optimum  beam  space  M=1 

A  Optimum  beam  space  M=2 

+  Optimum  beam  space  M=3 

X  Optimum  beam  space  M=4 

O  Optimum  beam  space  M~5 


Figure  4.  Array  gains  in  uncorrelated  receiver  noise 
plus  10  dB  interference  when  r/X  =  1/8 
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Optimum  beam 
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4- 

Optimum  beam 

space  M=10 

Figure  7.  Array  gains  in  isotropic  noise  plus  10  dB  interference 
plus  -10  dB  receiver  self  noise  at  r/A  =  1/8 
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Figure  8 


Figure  8.  Block  diagram  for  comparison  of  conventional,  array  space  and  beam  space  array  gain  estimates 


(SGI  NldO 


—  Conventional 
□  Optimum  array  space 

0  Optimum  beam  space  M=3 

A  Optimum  beam  space  M=6 

+  Optimum  beam  space  M=9 


Figure  9.  Array  gain  estimates  for  r/A 
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